Anomalous spin-waves and the commensurate-incommensurate magnetic phase 

transition in LiNiP04 
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Detailed spin-wave spectra of magneto-electric LiNiP04 have been measured by neutron scattering 
at low temperatures in the commensurate (C) antiferromagnetic (AF) phase with Tjv — 20.8 K. An 
anomalous low-energy mode is observed at the modulation vector of the incommensurate (IC) AF 
phase appearing above Tat. A linear spin- wave model based on Heisenberg exchange couplings and 
single ion anisotropics accounts for all the observed spin-wave dispersions and intensities. Along 
the b axis an unusually strong next-nearest-neighbor AF coupling competes with the dominant 
nearest-neighbor AF exchange interaction and causes the IC structure. 

PACS numbers: 75.30.Ds, 75.10.Jm, 75.50.Ee, 75.25.-l-z 
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Strongly correlated systems with electronic and mag- 
netic interactions have attracted much attention recently. 
Noted among them are the high-temperature supercon- 
ductors, the colossal-magneto-resistance materials, and 
the multiferroics [l[ . The latter group of materials possess 
magnetic and ferroelectric phases that may even co-exist 
over a range of temperatures 0]. For some, the correla- 
tions are manifested in the magneto-electric effect (ME) 
in which a magnetic field induces an electric polarization 
and an electric field induces a magnetic moment 0] . The 
lithium orthophosphatcs LiA/P04, with M = Mn, Fe, Co 
and Ni constitute a prototypical iso-structural group of 
antiferromagnets (AF) that are ME below their Neel tem- 
peratures fj, 0] . Although they have been investigated 
for long [J|, a general consensus about the microscopic 
interactions leading to their ME effect has not yet been 
established 0,01. 

Here we focus on LiNiP04 which compared to its 
lithium orthophosphate (LO) counterparts have more 
complex magnetic ordering properties and a different 
ternperature dependence of the linear ME coefficients 
[1, Q. LiNiP04 has an orthorhombic crystal structure 
with Prima (no. 62) symmetry, lattice parameters a = 
10.02 A, 6 = 5.86 A, c = 4.68 A, and four formula units 
per unit cell . The magnetic Ni^+ ions with spin S = 1 
are located on the 4(c) sites and form buckled planes per- 
pendicular to the crystallographic a axis. All LO com- 
pounds have low temperature zero-field commensurate 
(C) AF phases, but only LiNiP04 passes via a first order 
phase transition at T^v = 20.8 K into an incommensurate 
(IC) AF phase with ordering temperature T/c = 21.7 K. 



Above T/c short range correlations persist up to 40 K [8| . 
In a recent study we have determined the magnetic struc- 
tures and {H, T) phase diagram as function of tempera- 
ture and magnetic field up to 14.7 Tesla applied along the 
c axis . Here we also suggested a microscopic model ac- 
counting for the symmetry and temperature dependence 
of the measured ME coefficients 0]. The present inelas- 
tic neutron scattering study is aimed at determining the 
magnetic couplings in LiNiP04 to improve our knowledge 
of the magnetic phases. Although the C phase has the 
general ordering vector kc = (0,0,0) we find an anoma- 
lous shallow minimum in the magnon dispersion along 
the b direction at a modulation vector w (0, 0.12, 0) 
(in reciprocal lattice units). This is in the range of the 
IC ordering vector, k/c = (0, q, 0), where q is increasing 
with temperature between 0.07 < q < 0.155 Q. Analyzes 
of the spin-wave spectra show that the minimum at 
results from competing nearest-neighbor (NN) and next- 
nearest-neighbor (NNN) AF interactions along the b axis, 
and that the low temperature zero field AF structure is 
on the verge between C- and IC-order. 

A high quality single crystal with approximate dimen- 
sions 5x5x9 mm'^ and weight 0.4 g was used for this 
study. Inelastic neutron scattering measurements were 
performed on the RITA-II spectrometer (SINQ, PSI) by 
varying the incident energy at constant q to obtain dis- 
persions along (g,l,0) and (0,1 -I- q,0). The RITA-II uses 
a seven-blade analyzer and a final energy of 5 meV [l^l ■ 
Spinwaves along the three principal directions, ((7,1,0), 
(0,1 + q,0) and (0,1, g) were measured on the HBIA spec- 
trometer (HFIR, Oak Ridge) with 14.7 meV initial en- 
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FIG. 1; Measured (filled circles) spin- wave dispersion along three reciprocal directions, compared with a Holstein-Primakoff 
spin-wave calculation (full and dotted lines) as explained in the text. The insets show the corresponding measured and calculated 
intensities (in arbitrary units) as function of wavevector. The dispersions along (H,1,0) and (0,K,0) were measured at T = 2 K 
on RITA-II, while the dispersion along (0,1,L) was measured at T = 10 K on HBIA. 



ergy and were within uncertainties the same as those 
measured on RITA-II. More than an entire Brillouin zone 
in the foc-plane was also measured on the INS spectrom- 
eter (ILL), using a 47 channel multi-analyzer-detector 
(MAD) [ll[ set for a final energy of 30 meV. 

The spin-wave dispersions and integrated intensities 
measured at RITA-II and HBlA along the a-, h- and c- 
directions at low temperatures are shown in Fig. [T] The 
((7,1,0) and (0,1,9) dispersions behave as expected for an 
antiferromagnct with an energy gap at the zone center, 
but the dispersion along (0,H-(7,0) is unusual. It is almost 
a constant at low q with a shallow minimum at q k, 0.12, 
corresponding to the observed IC ordering vectors just 
below Tic @l- Generally, the (0,1 + 9,0) modes are lower 
in energy than along the two other principal directions, 
and they soften the same way as the 9 = mode when the 
temperature increases towards the C-IC transition tem- 



perature T/v. Fig. [2] shows the temperature dependence 
of the spin wave gap at Q = (0, 1,0). The inset displays 
spin-wave neutron scattering spectra at Q = (0, 1.055, 0), 
close to the g-value of the IC-modulation vector just 
above Tj^ . No spin- waves were observed in the IC phase. 

As a basis for our spin-wave model we adopt the AF 
ground state of the Ni^+ ions determined by Santoro 
et al. IJ] and Vaknin et al. [l3| with antiparallel spins 
pointing along the c-axis as shown in Fig.Oj- In fact, the 
spins are canted slightly away from the c axis con- 
sistent with a Dzyaloshinsky-Moria interaction, but this 
has negligible influence on the linear spin-wave model 
presented here. Accordingly, we use the following spin 
Hamiltonian for LiNiP04, with Heisenberg interactions 
Jij and standard single- ion anisotropy terms Da[S")'^ 
{a = a, 6, c): 



= ^JyS,,.S,+^i?„(5f)^ 
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FIG. 2: (a) Spin- wave energy gap at Q = (0, 1, 0) as function 
of temperature measured at RITA-II. The dashed line is a 
guide to the eye. The inset shows constant Q = (0, 1.055,0) 
scans at different temperatures with huj (in meV) along the 
a;-axis and measured intensity (in counts per minute) along y. 
(b) The magnetic unit cell of LiNiP04 showing only Ni^^ ions. 
The depicted spin configuration and exchange interactions are 
the ones used in the spin-wave model. 



In the 6c-plane the Ni-spins are coupled via super- 
exchange Ni^+-0-Ni^+ bonds, and our model includes 
NN interactions, Jha with a 3.81 A bond length, and two 
NNN, Jf, and Jc, with 4.71 A and 5.89 A bonds, respec- 
tively. For the couplings between the 6c-planes we con- 
sider NN interactions Jab and Jac, which have 5.40 A and 
5.50 A bond lengths (see Fig. Wp)- The inter-plane cou- 
plings arc mediated by P04-tctrahedra and may have a 
significant magnitude, e.g. as found in Li3Fe2(P04)3 [l^. 

Using linear spin-wave theory we have calculated 
the spin-wave dispersions and intensities by a Holstein- 
Primakoff approach similar to the one described in [l6| . 
The two spin-wave dispersions derived from Eq. ([1]) are 
given by (for details see [l4j] ) . 



hw = -{B± cy 



(2) 



where, 



3 



A = AS{Jbc + Jab) - 2S[Jb{l - cos(Q • rs)) + Jc(l - cos(Q ■ rg)) + Jac{2 - cos(Q • Yj) - cos(Q • rg))] + {S- 1/2) [D^ + Dt), 



B = {S-l/2)iDa-Db), C 



2 Jbc5'[cos(Q • ri) + cos(Q • r2)] 
I 



2Jq{,S'[cos(Q • rs) + cos(Q • r4)]. 



Here denote vectors between NN or NNN Ni-ions, and 
are: ri,2 = ^(b ± c), r3^4 = i(a ± b), = b, rg = c, 
r7,8 = ^(a±c). 

The model parameters were determined by a simulta- 
neous least squares refinement of the spin- wave disper- 
sions in all three principal directions [l7[. The best fit 
is shown by full and dotted lines in Fig. [TJ and the re- 
fined parameters are listed in Table [J The couplings Jac 
and Jc assure ferromagnetically ordered ac-planes that 
interact via J^c + Jab to NN and to NNN ac-planes 
along the ^-direction. All interactions Jf,c, Jab (NN) and 
Jb (NNN) are AF which leads to frustration promoting 
the IC-ordering. We note that the Jab and Jac couplings 
in the a-direction are somewhat weaker than Jbc and Jb- 
This is consistent with the findings in [ist where magnetic 
short range fluctuations extending into the paramagnetic 
regime and a critical exponent (3 = 0.12 of the C order 
parameter indicates a near-2D magnetic ordering of the 
&c-planes. The single ion anisotropics, Da and Db, are 
both positive indicating that a c axis magnetic moment is 
favored in the ground state, as observed experimentally. 
It should be noted that it is not possible from Eq. ^ to 
determine which of the two is larger. Da or Db- However, 
comparing the calculated intensities to the measured in- 
tensities (insets of Fig. [1]) determines unequivocally that 
Da < Db, as given in Table HI 

From the symmetry of the magnetic structure it can be 
shown that at least two non-degenerate dispersions are 
needed to account for the spin-wave in the 5-direction 
14| . In our model calculations these two non-degenerate 
branches is a result of the different anisotropics along a, 
b and c. However, at RITA-II and HBIA (Fig. [J) we ob- 
served only a single magnon dispersion, either because 
the second branch had too low intensity (see Fig. [T] in- 
sets), or because the two dispersions were indistinguish- 
able within instrumental resolution at the measured Q- 
values. Searching for the second branch and further con- 
firmation of our model, we used the INS spectrometer to 
collect data from a larger range of Q-values in the bc- 
plane. In Fig. Owe compare a 3D color map of the most 



TABLE I: The fitted spin coupling constants for LiNiP04. 
All units are in meV. By definition: Dc = meV. 



intense spin-wave dispersion measured at INS (a) to the 
results of our model calculation (b) using the interac- 
tion parameters of Table [J The experiment at INS gave 
clear evidence of two non-degenerate magnon branches, 
as seen in Fig. (4] Here two dispersions were observed at 
scattering vectors where the calculations predicted that 
the branches were well separated in energy and both had 
finite intensity. 

The exchange interactions established in this work may 
be used to explain the magnetic ground state. Maximiz- 
ing J(Q) = J2Jij^~^^ ^'' ^t zero temperature shows 
that IC-order is favorable in a simple model of layered 
magnetic systems with competing interactions if the ef- 
fective exchange interactions, Ji and J2 > 0, between 
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1.04(6) 0.670(9) -0.05(6) 0.30(6) -0.11(3) 0.339(2) 1.82(3) 
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FIG. 3: Spin-wave dispersions in the (0, K, L) plane corre- 
sponding to the most intense branch. The coloring is pro- 
portional to the observed spin-wave intensity (red for high 
intensities, blue for low), (a) Experimental data from INS at 
ILL. (b) Calculated data using the spin-wave model described 
in the text. The mesh around meV shows the deviation in 
energy between (a) and (b). 
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FIG. 4: (a) Constant Q = (0, 1, L) cuts from INS. At Q = (0,1,0.4) and (0,1,1.4) only a single branch is visible, but at 
Q = (0, 1, 1.8) and (0, 1, 2) both branches are observed, (b) Dispersions along (0, 1, L) measured at INS (blue and red circles) 
compared to the calculated spin- wave energies (full and dotted lines). The inset shows the calculated intensity of the dispersions. 
Two branches are observed at Q where the model predicts spin-waves with finite intensity that are well separated in energy, (c) 
Dispersions along (0, K, 2) measured at INS and compared to the calculated spin-wave energies. The inset shows the calculated 
intensity of the branches. Along {0,K,0) (Fig. [T]^) the intensity of the second branch is always zero, but here, along (0,-?^, 2), 
there are Q- values where both branches have intensity and can be observed. 



NN and NNN ferromagnetic layers fulfill that | Ji | < 4 J2 
[lit . If tliis condition applies, the magnetic ordering vec- 
tor Q is determined by cos{Qd) = — J1/4J2, where d is 
the inter-layer distance between adjacent ferromagnetic 
layers. Since LiNiP04 has ferromagnetic layers perpen- 
dicular to the ^-direction, we can test this condition. Here 
Ji = 2Jbc + 2Ja6 and J2 = Jt, while d = 6/2. Using 
the exchange parameters in Table H] we find: Ji = 2.7(2) 
meV and 4J2 = 2.68(4) meV, for the competing inter- 
actions. Witliin the uncertainties the magnetic ground 
state could therefore be either C or IC, and the system 
is close to an instability. A mean-field model using inter- 
action parameters obtained from the spin-wave spectra 
predict a behavior qualitatively similar to that observed 
experimentally [l^. Here it is found that the magnetic 
structure is IC just below Tic, but that the extra lock- 
in anisotropy energy gained by the C structure becomes 
important at slightly lower temperatures and as result 
changes q from k/c to 0. 
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